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ABSTRACT 

The relationship between the black hole mass and velocity dispersion indicated with [O III] 
line width is investigated for a sample of 87 flat-spectrum radio quasars (FSRQs) selected 
from SDSS DR3 quasar catalogue. We found the il/bh — f [oiii] relation is deviated from 
Tremaine et al. relation for nearby inactive galaxies, with a larger black hole mass at given 
velocity dispersion. There is no strong evidence of cosmology evolution in A/bh — cr[oiii] 
relation up to z ~ 0.8. A significant correlation between the [O III] luminosity and Broad 
Line Region (BLR) luminosity is found. When transferring the [O III] luminosity to Narrow 
Line Region (NLR) luminosity, the BLR luminosity is, on average, larger than NLR one by 
about one order of magnitude. We found a strong correlation between the synchrotron peak 
luminosity and NLR luminosity, which implies a tight relation between the jet physics and 
accretion process. 
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1 INTRODUCTION 

The evolution of black holes has been shown to be closely coupled 
to that of their host galaxies for normal galaxies, mainly through 
the tight correlation between the central black hole mass and the 
bulge mass and luminosity (Kormendy & Richstone 1995; Magor- 
rian et al. 1998), and tighter one between the central black hole 
mass and stellar velocity dispersion a, in the galactic bulge (Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000a). The latter relation 
has been established by Tremaine et al. (2002) for a sample of 31 
nearby inactive galaxies as 
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For some AGNs with available bulge velocity dispersion and the 
reverberation mapping black hole mass, Gebhardt et al. (2000b) 
and Ferrarese et al. (2001) found that these AGNs also follow the 
Mhh — cr* relation founded in the nearby inactive galaxies. With 
the aim to constrain the nature and evolution of AGNs and the ad- 
vantage of relative higher redshift of AGNs, the i\fbh ~ cr» rela- 
tion and its evolution have been extensively explored for different 
AGN populations, such as radio-quiet, radio-loud AGNs, narrow 
line Seyfert 1 galaxies and young radio galaxies et al. (e.g. Shields 
et al. 2003; Bian & Zhao 2004; Bonning et al. 2005; Liu & Jiang 
2006; Salviander et al. 2007; Bian et al. 2008; Shen et al. 2008; Wu 
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2009a), however it is still in large debates whether AGNs generally 
follow the Afbh — o"» relation. 

While the derivation of A/bh from AGNs broad line widths 
and continuum (and/or broad line) luminosity is now well estab- 
lished (e.g. Kaspi et al. 2000; Greene & Ho 2005; Kong et al. 2006), 
the measurements of cr, are generally difficult for quasars, due to 
the faintness of the host galaxy and the relative brightness of active 
nucleus. When stellar velocities dispersion cannot be measured, the 
line widths of the narrow emission lines (e.g. [O III], [S II]) usu- 
ally can be used to trace a, (Nelson & Whittle 1996; Greene & Ho 
2005). Nelson & Whittle (1996) compared the bulge magnitudes, 
[O III] A5007 line widths, and stellar velocity dispersions in Seyfert 
galaxies, adopting (T[oiii] = FWHM([OIII])/2.35, and found, on 
average, good agreement between (J[oni] and a*, although (T[oiii] 
shows more scatter than cr» on a Faber-Jackson plot. This implies 
that the kinematics of narrow-line region (NLR) gas is largely dom- 
inated by the bulge gravitational potential, therefore, can be effec- 
tively used as a substitute for a, of galaxy bulges. Using [O III] 
line width as surrogates for cr,, Nelson (2000) claimed that the 
Afbh — o"» relation for normal galaxies and active galactic nuclei is 
preserved. Albeit some defect of [O III] line profile (e.g. asymme- 
try and non-Gaussian) and other surrogates proposed (e.g. [S II], 
Greene & Ho 2005), the width of [O III] is still most commonly 
used as surrogate for ct« in AGNs (e.g. Bonning et al. 2005; Bian 
et al. 2008). 

The formation of radio jets is still a unresolved issue in AGNs 
research. While radio-quiet AGNs roughly follow the Afbh ^ 
relation of normal galaxies, radio-loud ones deviate, in the sense 
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that radio-loud objects have, on average, larger black hole mass 
than radio-quiet objects for a given velocity dispersion (e.g. Bian 
&. Zhao 2004; Bonning et al. 2005; Liu & Jiang 2006; Salviander 
et al. 2007; Bian et al. 2008; Shen et al. 2008). Among the popula- 
tion of radio-loud AGNs, FSRQs represent an extreme class, which 
are generally characterized by strong and rapid variability, high po- 
larization, and apparent superluminal motion. These extreme prop- 
erties are generally interpreted as a consequence of non-thermal 
emission from a relativistic jet oriented close to the line of sight. 
In this letter, we investigate the relationship between the black hole 
mass and o-[oiii] for a sample of 87 SDSS FSRQs based on the 
spectral analysis of a larger FSRQs sample in Chen et al. (2009). 
The cosmological parameters Ho = 70kms~^ Mpc~^, f2ni=0.3, 
Q.A = 0.7 are used throughout the paper, and the spectral index a is 
defined as fi, oc i/^" with being the flux density at frequency v. 



2 SAMPLE AND DATA ANALYSIS 

The parent sample of this work is 185 FSRQs in Chen et al. (2009). 
It was constructed through cross-correlating the Shen et al. (2006) 
SDSS DR3 X-ray quasar sample (3366 sources, see Shen et al. 
2006 for details) with Faint Images of the Radio Sky at Twenty- 
Centimeters 1.4 GHz radio catalogue (FIRST, Becker, White & 
Helfand 1995) and the Green Bank 6 cm radio survey at 4.85 GHz 
radio catalogue (GB6, Gregory et al. 1996). The sample of 185 FS- 
RQs was constructed from conventional definition of FSRQs with 
a spectral index between 1.4 and 4.85 GHz a < 0.5. In this work, 
we select 87 FSRQs with redshift z < 0.83 out of 185 FSRQs. The 
redshift restriction is selected in order to measure [O III] 5007A 
line from SDSS spectra as well as H/3. 

The SDSS spectra cover the wavelength range from 3800 to 
9200 A with a resolution of about 1800 - 2000 (see Schneider et 
al. 2005 for details). The spectral analysis is briefly described in 
this work, and more details can be found in Chen et al. (2009). The 
SDSS spectra were firstly corrected for the Galactic extinction us- 
ing the reddening map of Schlegel, Finkbeiner & Davis (1998) and 
then shifted to their rest wavelength. We choose those wavelength 
ranges as pseudo-continua, which are not affected by prominent 
emission lines, and then decompose the spectra into the follow- 
ing three components: (1) A power-law continuum to describe the 
emission from the active nucleus. (2) An Fe II template adopting 
the UV Fe II template from Vestergaard & Wilkes (2001), and op- 
tical one from Veron-Cetty et al. (2004). (3) A Balmer continuum 
generated in the same way as Dietrich et al. (2002). 

The modeling of above three components is performed by 
minimizing the in the fitting process. The final multicomponent 
fit is then subtracted from the observed spectrum. The broad emis- 
sion lines were measured from the continuum subtracted spectra. 
For the redshift range of our sources, we focused on several promi- 
nent emission lines, i.e. Hq, H/3 and Mg II. Generally, two gaussian 
components were adopted to fit each of these lines, indicating the 
broad and narrow line components, respectively. The blended nar- 
row lines, e.g. [O III] AA4959, 5007A and [He II] A4686A blend- 
ing with H/3, and [S II] AA6716, 6730A, [N II] AA6548, 6583A 
and [O I]A6300A blending with Ha, were included as one gaus- 
sian component for each line at the fixed line wavelength. The 
minimization method was used in fits. The line width FWHM, line 
flux of broad Ha, H/3, Mg II and narrow [O III]5007A lines were 
obtained from the final fits for our sample. 

There are various empirical relations between the radius of 
broad line region (BLR) and the continuum (and/or broad line) lu- 



minosity, which can be used to calculate the black hole mass in 
combination with the line width FWHM of broad emission lines. 
Since the continuum luminosity of FSRQs are usually contam- 
inated by the non-thermal jet emission, we use the Vestergaard 
& Peterson (2006) method to calculate Mbh, which utilizes the 
FWHM and luminosity of broad H/3: 

In this work, we also calculate the BLR luminosity Lblr fol- 
lowing Celotti, Padovani & Ghisellini (1997) by scaling the strong 
broad emission lines Ha, H/3 and Mg II to the quasar template 
spectrum of Francis et al. (1991), in which Lya is used as a ref- 
erence of 100. By adding the contribution of Ha with a value of 
77, the total relative BLR flux is 555.77, of which H/3 22 and Mg 
II 34 (Celotti et al. 1997; Francis et al. 1991). 

As shown in Chen et al. (2009), the spectral energy distri- 
bution (SEDs) were constructed for each source using the radio 
FIRST 1.4 GHz and GB6 4.85 GHz data, the optical data selected 
as the line-free spectral region from SDSS spectra, and IkeV X- 
ray data compiled in Shen et al. (2006) from ROSAT All Sky Sur- 
vey. When available, the 2MASS IR (J, H, and Ks) (Skrutskie 
et al. 2006) and the Far- and near-UV GALEX data (Martin et 
al. 2005) are also added. The synchrotron peak frequency and the 
corresponding peak luminosity were thus obtained for each source 
through fitting the SED with a third-degree polynomial following 
Fossati et al. (1998). The thermal emission from accretion disk and 
host galaxy is estimated in Chen et al. (2009). The accretion disk 
emission is calculated assuming a steady geometrically thin, opti- 
cally thick accretion disk around a Schwarzschild black hole and 
using estimated black hole mass and bolometric luminosity. The 
host galaxy spectra is estimated using the elliptical galaxy template 
of Mannucci et al. (2001), in combination with the bulge abso- 
lute luminosity in R-band estimated from Mbh — Mr relation of 
McLure et al. (2004). The total contribution of accretion disk and 
host galaxy thermal emission are estimated by calculating the frac- 
tion of the thermal emission to the SED data at SDSS optical and 
GALEX UV region. A marginal value of 50% at most of SED 
wavebands is used to divide the FSRQs into thermal-dominated 
(> 50%) and nonthermal-dominated (< 50%). Using this crite- 
rion, 36 of 87 FSRQs in this work are non-thermal dominant FS- 
RQs, and 51 are thermal-dominated ones. The [O III] measure- 
ments are finally performed in 81 FSRQs, and [O III] line is rather 
weak in remaining 6 FSRQs, which therfore can not be properly 
measured. Among 81 [O III] available FSRQs, 32 sources are non- 
thermal dominant FSRQs, and 49 are thermal-dominated FSRQs. 



3 RESULTS AND DISCUSSION 

3.1 Mbh - cr[oiii] relation 

The relationship between the black hole mass Mbh and velocity 
dispersion (j[oiii] = FWHM([OIII]5007A)/2.35 for all 81 FS- 
RQs is shown in Fig. 1, in which the Tremaine et al. relation of 
equation (1) is also plotted as dashed line. From the figure, the sig- 
nificant scatter and deviation from Tremaine et al. relation is ap- 
parently seen, with majority of sources lying above the relation, 
which implies a higher black hole mass for our FSRQs than nearby 
inactive galaxies at given velocity dispersion. This result is actu- 
ally consistent with Liu & Jiang (2006) and Bian et al. (2008) for 
general radio-loud AGNs/quasars. Both works claimed that radio- 
quiet AGNs evenly follow the Tremaine et al. Mbh — cr* relation. 
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Figure 1. The velocity dispersion o"[oiii] versus black hole mass. The 
dashed line is the M|-,h — relation of Tremaine et al. (2002) for 3 1 nearby 
inactive galaxies. The solid circles are non-thermal dominant FSRQs, while 
the open circles are thermal-dominant FSRQs. 
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Figure 2. The luminosity at 5100A versus broad H/3 luminosity. The sym- 
bols are same as in Fig. 1 . The dashed line is the OLS bisector linear fit to 
Kaspi et al. (2000) radio quiet AGNs (see Liu et al. 2006). 



however, radio-loud ones deviate with the same trend as ours. Al- 
ternatively, radio-loud quasars are found to have relatively smaller 
[O III] line width than radio quiet quasars at given black hole mass 
or bulge luminosity (e.g. Shields et al. 2003; Bonning et al. 2005). 
However, the reason is still unknown. 

The relation between the BLR radius and H/3 luminosity, and 
the corresponding empirical relation to estimate black hole mass in 
equation (2), are basically scaled from the AGNs sample with re- 
verberation mapping black hole mass, of which most of sources 
are radio-quiet. The main advantage of equation (2) in estimat- 
ing black hole mass for FSRQs is that the contamination of non- 
thermal jet emission (likely Doppler boosted) in continuum lumi- 
nosity can be avoided. However, it is not proved yet whether it is 
applicable to radio-loud AGNs. As an evaluation of jet emission 
contamination, Liu et al. (2006) shown that the optical continuum 
luminosity of FSRQs usually exceed the value estimated from the 
Lhs - L relation L = = 0.843 x 10^ L?,?® fitted for 

5100A 5100A 

Kaspi et al. (2000) radio quiet AGNs using ordinary least-square 
(OLS) linear fit method (see Fig. 4 in Liu et al. 2006). This result 
is consistent with the expectation that the continuum emission of 
FSRQs are contaminated by the (even dominant) non-thermal jet 
emission. As a check for our FSRQs sample, we plot the relation- 
ship between Lh/j and L^^^ ^ in Fig. 2. We found that the ther- 
mal FSRQs well follow the tup — L^^^ ^ relation for radio quiet 
AGNs of Liu et al. (2006), supporting t^ie dominance of thermal 
emission in these 51 FSRQs. However, large scatter exist for 36 
nonthermal FSRQs. While most sources follow the relation, about 
10 of 36 sources deviate from the relation, with larger luminosity 
at 5100A than expectations from relation, which is most likely due 
to the contribution of non-thermal jet emission. 

The consistence of Lhs ~ L - relation of the thermal- 

5100A 

dominant FSRQs with that of radio quiet AGNs implies that the 
thermal continuum emission and the broad line emission follow 
the similar relation. Therefore, there will be likely little problem 
in estimating the BLR radius using empirical Lh/j — i?BLR re- 
lation, unless the BLR kinematics (e.g. BLR radius) in our FS- 
RQs is largely different from that of radio quiet AGNs. However, 
it is rather difficult to move our sources downward to follow the 
Tremaine et al. relation solely by changing BLR radius, let alone 



to reduce the large scatter in Fig. 1 . For radio-loud quasars, espe- 
cially FSRQs, the disk-like BLR geometry (Wills & Browne 1986) 
can cause the underestimate of black hole mass (Lacy et al. 2001; 
McLure & Dunlop 2001). However, this effect will even make sit- 
uation more sever, i.e. our FSRQs will deviate towards larger black 
hole mass from the Tremaine et al. relation. The likely explanation 
of smaller [O III] line width in radio-loud quasars can be a different 
[O III] kinematics or geometry in radio-loud quasars, especially for 
our FSRQs. However, the reason is unclear. The clues, if any, may 
be from the recent finding that the [O III] luminosity is not emit- 
ted isotropically and that there is significant extinction towards or 
within the narrow line region in a subset of Seyfert 2 galaxies, in the 
way that the observed [O III] luminosity are systematically smaller 
for obscured Seyferts by comparing [O III] with [O IV] 25.9 /^m for 
a unbiased Seyfert galaxies sample (Diamond-Stanic et al. 2009; 
see also Jackson & Browne 1990; Hass et al. 2005; Melendez et 
al. 2008; Zhang et al. 2008). More relevant to our FSRQs as radio 
powerful sources, the [O III] emission from quasars was compared 
with radio galaxies by Jackson & Browne (1990), and they found 
that the [O III] emission of quasars is much stronger than that of ra- 
dio galaxies, which was interpreted as part of the [O III] emission 
being obscured by the torus in radio galaxies. The similar conclu- 
sion was reached by Haas et al. (2005) by comparing [O III] with 
[O IV] 25.9 /im for a sample of seven 3CR FR II radio galaxies and 
seven 3CR quasars. It is well established that the NLR is stratified 
with high-density and high ionization gas at close to the continuum 
source whereas low density and low-ionization gas is in the outer 
part of the NLR (Nagao et al. 2003; Riffel et al. 2006). Zhang et al. 
(2008) claimed that a significant fraction of the [O III] emission, at 
least for their Seyfert galaxies, may arise from the inner dense part 
of the NLR, which however can be covered thus obscured by the 
torus with its inner edge on scales of parsecs, and its extent likely 
on scales of several tens of parsecs (e.g., Schmitt et al. 2003, Jaffe et 
al. 2004). Due to the complex of [O III] line profile (Bonning et al. 
2005), nevertheless, it needs further investigation on the Mbh — o"« 
relation for our FSRQs using other narrow lines, e.g. [S II] (Greene 
& Ho 2005). However, the redshift range will be much limited. 

Although the detailed reason is unclear, the deviation of our 
FSRQs in Fig. 1 could be related with several recent findings on 
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the bulge - black hole relations (e.g., Aller & Richstone 2007; Hop- 
kins et al. 2007; Lauer et al. 2007). The black hole masses pre- 
dicted from the Mbh — o relation are found to be in conflict with 
those from the Mbh — ^buigo relation for high-luminosity galax- 
ies (e.g. brightest cluster galaxies: BCGs), with the former relation 
predicting a larger black hole mass (Lauer et al. 2007). While this 
may be explained by the slow increase in o with Lbuigc and the 
more rapid increase in effective radii with Lbuigc seen in BCGs 
as compared to less luminous galaxies, the authors argued that the 
Mbh — ibuige relationship is a plausible description for galaxies of 
high luminosity. From the major galaxy merger simulations, Hop- 
kins et al. (2007) found the so-called black hole mass fundamental 
plane in the form of Afbh oc m^'^^^^-^"^ a^-'^^°-^ , similar to re- 
lations found observationally, where A/, is the stellar mass. They 
also claimed that this fundamental plane is better than any single- 
variate predictor of black hole mass, e.g. the Mbh — o" relation and 
the Mbh — M, relation. Moreover, the bulge gravitational binding 
energy may play a dominant role, in that A/bh oc iSg '' as found 
by Aller & Richstone (2007), which is as strong a predictor of Afbh 
as the velocity dispersion a, for the elliptical galaxies. In view of 
these findings, it is possible that our FSRQs, especially those with 
high black hole mass (e.g. A/bh > 3 x lO^M©), may have larger 
bulge mass at a given cr. Alternatively, if these objects lie in the cen- 
ter of galaxies that continue to accrete the surrounding ICM after 
the initial black hole growth, as in so-called radio-mode feedback 
models (e.g., Croton et al., 2006), then perhaps the black hole can 
grow to a higher mass than the mass when it originally settled on 
the Mbh — o" relationship. As shown in Croton et al. (2006) (see 
their Fig. 3.), the growth of black holes is dominated by the 'quasar 
mode' at high redshift and falls off sharply at z; < 2, in contrast, 
the 'radio mode' becomes important at low redshifts where it sup- 
presses cooling flows. 

The deviation of A/bh from Tremaine et al. relation can be cal- 
culated as AA'fbh = Mbh - Mbh(cr[oiii]), of which Mbh(o-[oiii]) 
is the predicted black hole mass from equation (1) using [O 
III] velocity dispersion. We perform the correlation analysis be- 
tween AMbh and other parameters, i.e. redshift, Eddington ratio 
iboi/iEdd, 5 GHz lumiuoslty I/sGHz, and synchrotron peak lumi- 
nosity Lpeak, and only found a moderately significant correlation 
between AA/bh and redshift with Spearman correlation coefficient 
r = 0.296 at confidence level ~ 99.3%. However, this correlation 
may likely be caused by the common dependence of other param- 
eters, e.g. I/H/3. Indeed, the partial Spearman correlation analysis 
shown that the strong correlation no longer exist independent of 
Lh/3- Therefore, there is no strong evidence of cosmology evolu- 
tion of Afbh — o"[oiii] relation up to z = 0.83. This is in contrast 
to the results of several studies (e.g. Woo et al. 2006; Shields et al. 
2006) that have found apparent positive correlations between the 
Mbh — cr* relation and redshift. However, Shields et al. (2003) and 
Salviander et al. (2006), found no evidence for redshift evolution 
up to 2 ~ 3, using [O III] line widths as surrogates for a* . 

3.2 Synchrotron peak luminosity, broad and narrow line 
region luminosity 

The broad and narrow emission lines are believed to be produced 
from photonionization process, although the shock-excitation from 
jet-ISM interaction can not be ignored for narrow emission lines 
in some AGNs (e.g. Dopita & Sutherland 1995). Therefore, the in- 
timate relation is expected between broad and narrow line lumi- 
nosities. Indeed, a significant correlation between the [O III] lumi- 
nosity and BLR one is found with Spearman correlation coefficient 
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Figure 3. The BLR luminosity versus NLR luminosity. The symbols are 
same as in Fig. 1. The solid line is the OLS bisector linear fit for all 81 
FSRQs. 

r = 0.875 at confidence level of S> 99.99%. Using partial Spear- 
man correlation analysis to exclude the common dependence of 
both luminosities on redshift, the correlation is still significant with 
correlation coefficient r — 0.664 at confidence level of S> 99.99%. 

In this work, we tentatively transfer the [O III] 5007A lumi- 
nosity to NLR luminosity, using relation of Lnlr = 3 x (3 x 
L[oii] + 1-5 X -L[oiii]) and assuming L[oiii] = 4 x Ljqii] (Rawl- 
ings & Saunders 1991). The relationship between the luminosity 
of BLR and NLR is presented in Fig. 3. The ordinary least-square 
(OLS) bisector linear fit to Lblr - Lnlr relation gives: 

log Lnlr = (0.89 ± 0.05) log Lblr + (3.91 ± 2.10) (3) 

The mean value of log Lnlr is (log Lnlr) = 43.52 ± 0.66, and 
for log Lblr is (log Lblr) = 44.53 ± 0.75. The BLR luminos- 
ity is, on average, larger than NLR luminosity by about one order 
of magnitude with (log (Lblr/Lnlr)) = 1-01 ± 0.37. This re- 
sult means that the covering factor of NLR is about one tenth of 
that of BLR. Moreover, it is consistent within a factor of three with 
the ratio of BLR to NLR luminosities (~ 25) calculated by us- 
ing the relative flux in the composite quasar spectrum of Francis et 
al. (1991), which shows that our FSRQs spectra generally follow 
the composite spectrum although they are radio powerful sources. 
In addition, the scaling between the BLR and NLR luminosities 
enables us to statistically estimate the bolometric luminosity from 
narrow emission lines, which can be more readily observed, e.g. 
for type 2 AGNs. If using Lboi ~ 30Lblr newly calibrated in Xu, 
Cao & Wu (2009), statistically, the bolometric luminosity can be 
estimated as Lboi = 300Lnlr, i.e. Lboi ~ 2025 L[oni]> imply- 
ing a covering factor of ~ 0.003 for narrow line region, similar to 
those of Rawlings & Saunders (1991) (~ 0.01) and Willott et al. 
(1999) (~ 0.003). More specifically, the calibration between the 
bolometric luminosity and NLR luminosity can be given through 
OLS bisector linear fit: 

log Lboi = (1.12 ± 0.06) log Lnlr - (3.92 ± 2.59) (4) 

Our bolometric correction is consistent with Heckman et al. (2004) 
(Lboi/L([OIII]) ~ 3500 for Seyfert galaxies with a scatter of 
approximately 0.38 dex) and Willott et al. (1999) (Lboi ~ 5 x 
10'^L[oii] W). The consistent results are also recently found by Wu 
(2009b) that Lboi oc Lfofj^" "* for a sample of Seyfert 1 galax- 
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Figure 4. The upper panel: The synchrotron peak luminosity versus NLR 
luminosity. The symbols are same as in Fig. 1. The solid line is the OLS 
bisector linear fit for 32 non-thermal dominant FSRQs. The lower panel: 5 
GHz luminosity versus NLR luminosity. The symbols are same as in Fig. 1. 



ies and radio-quiet quasars using Lboi = 9ALa(5100A), and the 
mean bolometric correction Lboi/i([OIII]) ~ 3400, 2000 for 23 
radio-quiet quasars and 20 Seyfert 1 galaxies, respectively. 

The relationship between the NLR luminosity and synchrotron 
peak luminosity is given for all 81 FSRQs in Fig. 4. We found a 
strong correlation between two parameters with a Spearman corre- 
lation coefficient r = 0.735 at confidence level ^ 99.99%. The 
partial correlation analysis gives a significant correlation with cor- 
relation coefficient r = 0.408 at confidence level ~ 99.98% inde- 
pendent of redshift. When only 32 non-thermal dominant FSRQs 
are considered, the original significant correlation (r = 0.637 at 
confidence level ~ 99.96%) still remains (r = 0.406 at confidence 
level ^ 98%). For these 32 non-thermal dominant FSRQs, the OLS 
bisector linear fit on i^L^pcak - -t/NLR relation gives: 

log i^i.pcak = (0.80 ± 0.10) log Lnlr + (10.46 ± 4.32) (5) 

which remains same when fitting on all 81 FSRQs. The strong cor- 
relation is also found between the NLR and 5 GHz luminosity (see 
Fig. 4), however, it is less significant (r = 0.463 at confidence level 
^ 99.99%) than that between the NLR luminosity and synchrotron 
peak luminosity. 



The advantages of using the synchrotron peak luminosity is 
that the most of synchrotron emission are radiated at synchrotron 
peak frequency, at which the luminosity can be a good indicator of 
synchrotron emission. Since the synchrotron peak frequency varies 
from source to source, the luminosity at fixed waveband is actually 
from the different portion of source SED. This can be strength- 
ened from Fig. 4, in which the large scatter in vLvzg — Lnlr is 
significantly reduced in i^Li/pcak — Lnlr, and with the less sig- 
nificance of the former correlation. This not only implies a tight 
relation between jet physics and accretion process, but also claims 
that the synchrotron peak luminosity can be a better indicator of 
jet emission than 5 GHz one. However, the possibility that part, if 
not all, of the i/L^poak — Lnlr correlation can be caused from 
shock-excitation for NLR by the jet-ISM interactions (e.g. Dopita 
& Sutherland 1995; Bicknell et al. 1997), can not be completely 
excluded. By comparing the [O IV] 25.9 /im and [O III] luminosi- 
ties for the combined sample of radio loud and radio quiet Seyfert 
galaxies, Melendez et al. (2008) found that radio loud sources ex- 
hibit higher emission-line luminosities than those of radio quiet 
ones. The authors claimed that this result could be explained by 
a proposed bow shock model where part of the NLR emission is 
being powered by radio-emitting jets, present in radio-loud AGNs. 
Moreover, the jet-ISM interaction are invoked to explain the highly 
broadened [O III] line profiles in young radio galaxies (CSS/GPS 
sources) and some luminous linear radio sources (e.g., Gelderman 
& Whittle 1994; Nelson & Whittle 1996). 

Although it can be a good indicator of the jet emission, the 
defect of the synchrotron peak luminosity lies in the contamina- 
tion from beaming effect, which precludes it to well indicate the 
intrinsic source power. Only when the Doppler boosting is known 
for each source, the intrinsic source power can be obtained, which 
however can not be performed at present stage. 
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